Tunneling density of states as a function of thickness in superconductor/ strong 

ferromagnet bilayers 
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We have made an experimental study of the tunneling density of states (DOS) in strong ferro- 
magnetic thin films (CoFe) in proximity with a thick superconducting film (Nb) as a function of 
dF, the ferromagnetic thickness. Remarkably, we find that as dp increases, the superconducting 
DOS exhibits a scaling behavior in which the deviations from the normal-state conductance have a 
universal shape that decreases exponentially in amplitude with characteristic length d* w 0.4 nm. 
We do not see oscillations in the DOS as a function of d_F, as expected from predictions based on 
the Usadel equations, although an oscillation in T c (dp) has been seen in the same materials. 



One of the incompletely solved problems in conven- 
tional (noncuprate) superconductivity is the interaction 
between superconductivity and magnetism. This issue 
arises most prominently in the context of the so-called 
magnet superconductors (e.g., CeCoIn 5 ) and in the su- 
perconductor/ferromagnet (SF) proximity effect. One 
striking effect expected for a superconductor in the pres- 
ence of an exchange field is the existence of spatial mod- 
ulations of the superconducting pair wave function (see 
for instance Ref. [jj]). These oscillations occur in a new 
superconducting state where the center of mass of pairs 
acquires a non-zero momentum. This state was predicted 
40 years ago and is known as the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state HQ. In the case 01 the SF 
proximity effect, this oscillating pair wave function is ex- 
pected to exponentially decay into the F layer. These 
oscillations, in turn, are predicted to lead to oscillations 
in the critical temperature of SF bilayers inversions 
of the DOS |EI3> an d changes in the sign of the Joseph- 
son coupling in SFS sandwiches || (creating a so called 
7r-junction), as the F layer thickness, dp, is varied. Even 
more exotic predictions include possible odd-frequency 
triplet superconductivity 0- 

Indeed, a vast theoretical literature now exists re- 
garding the SF proximity effect, mostly concerning sys- 
tems with a uniform magnetization in the F layer in the 
dirty quasiclassical limit (i.e. where all the characteristic 
lengths are larger than both A^, the Fermi wavelength, 
and £, the mean free path). In this situation, the super- 
conducting properties can be calculated using the Usadel 
equation [l(J. Still, explicit predictions may differ de- 
pending on the importance of scattering processes Q, [~ 
or on the boundary conditions, which can be resistive 
or magnetically active |12j . 

Experimentally, phase sensitive measurements in some 
SFS structures clearly demonstrate the existence of ir- 
junctions 0,0,0]. On the other hand, critical temper- 
ature, T c , measurements in SF structures have shown a 
variety of behaviors as a function of dp . Everything from 



monotonic dependence to step-like featu res, small dip s 
and oscillations have been reported pi IwL []J\ . IfSl fl9t l20j . 

It has been pointed out that important parameters such 
as the size of the exchange field, E ex , and the bound- 
ary resistance, jb, can evolve naturally as a function of 
dp |23. What makes these results particulary hard to 
interpret is that all of these changes typically take place 
within a few nanometers, just like the expected oscilla- 
tion of the superconducting wave function inside the F 
material. 

In short, the situation is complicated both theoreti- 
cally and experimentally, and no clear, comprehensive 
understanding has emerged. This suggests that some 
new experimental approaches that probe superconduc- 
tivity directly inside the F material might be helpful. 

In this paper we present density of state spectroscopy 
studies using tunneling junction located on the F side 
of Nb/Coo.6Feo.4 bilayers. CoFe is a strong ferromag- 
net with a Curie temperature of approximately 1 100 K 
widely used in magnetic tunnel junction devices due to 
the high quality of the interface it makes with aluminum 
oxide, which is now the standard choice for tunneling 
barriers. The critical temperatures of similarly deposited 
Nb/CoFe bilayers, but with a thinner Nb layer (18 nm), 
show an oscillatory behavior as a function of dp: a slight 
dip is noticeable before saturation at large dp- A quan- 
titative analysis of these data based on the Usadel equa- 
tions appears elsewhere [i^- Based on these data, we 
expect that any inversion in the DOS that may arise in 
our samples would occur at the same thickness as the dip 
in T c , which is between 1 and 2 nm. Thus, we performed 
tunneling spectroscopy on samples with thicknesses rang- 
ing from to 4.5 nm in increments of 0.5 nm. 

Our results are remarkable in that we find that the 
deviations of the density of states from that of the normal 
state exhibit a precise scaling in amplitude over many 
orders of magnitude as dp is increased. No oscillations 
in the sign of these deviations are observed. 

Our junctions are deposited and patterned entirely in 
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FIG. 1: Raw conductance data taken at 0.5 K in a 
Nb/CoFe/AlOx/Al/CoFe junction at zero field and above the 
Nb critical field. Here cLf = 2.5 nm. Left inset: sample top 
view. Right inset: junction cross section. The black areas 
represent Al oxide layers. 



site in a DC magnetron sputtering system, described here 
|22| . The full geometry of the tunneling structures is pic- 
tured in the inset of Fig. ^ A thin (3 nm) CoFe layer 
which suppresses superconductivity in the Al electrode 
is deposited first — we see no superconducting transi- 
tion in the Al above 0.3 K — then an atomic oxygen 
source is used to fully oxidize the barrier to a thickness 
of ss 2 nm. The stencil mask is then changed and without 
breaking vacuum a CoFe layer and then a Nb layer are 
immediately deposited, followed by a very thin Al cap to 
protect the Nb layer from oxidation. In order to establish 
a robust superconducting state in the S layer, we deposit 
approximately 50 nm of Nb, which is thick enough that 
there should be only small changes in T c (dp) |30j. Typi- 
cal junction resistances are roughly 200 fi, with a trend 
to higher junction resistances with increased CoFe thick- 
ness. 

Measurements shown here are taken at 0.5 K using 
standard lock-in techniques to measure the differential 
resistance; variations as small as few parts in 10,000 could 
be distinguished. This resolution is crucial since when dp 
exceeds 1 nm the superconducting signal becomes very 
weak, reaching less than a part per thousand at dp — 
3 nm. 

Figure shows typical conductance spectra taken at 
0.5 K in zero-field and with a perpendicular field of 2.5 T, 
just above the critical field, H C 2. The zero- field curve 
clearly reveals the superconducting density of states, but 
it is superposed on top of a zero-bias anomaly. The high 
field data show the zero-bias anomaly on its own. This 
anomaly is a V-shaped curve, centered at zero volts, with 
the conductance at 1 mV typically two percent larger 
than the zero-bias conductance. The main experimen- 
tal difficulty lies in extracting the superconducting DOS 
from this zero-bias anomaly. To do this, we divide the 
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FIG. 2: Normalized conductances taken at 0.5 K for various 
CoFe thickness indicated inside. From the bottom to the top 
plot the vertical scale is successively amplified. The curves 
are shifted for clarity. The dashed line shows that the peaks 
maximum remains unchanged from 0.5 to 3 nm. 



zero-field curve by the in-field curve. This straightfor- 
ward procedure alone was sufficient for the thicker bar- 
riers, because no change in the normal conductance was 
seen between perpendicular and parallel applied fields or 
as a function of field strength. However for samples with 
thin CoFe layers (less than 1.5 nm), an additional zero- 
bias resistance peak increasing with the applied field was 
seen above H C 2, similar to previous studies [23(. For these 
samples, the zero-field normal-state background could be 
calculated by studying the field-dependence of this fea- 
ture and extrapolating it to zero-field. 

Figure [3 shows the resulting curves for all measured 
thicknesses. Note that in the top panel the scale is am- 
plified about a thousand times. Beginning at the bottom 
of the figure, we see that for dp — a clean DOS curve for 
Nb is obtained, as expected. A BCS fit to this curve gives 
1.3 meV for the gap energy. When the thinnest (0.5 nm) 
F layer is added to the superconductor, the tunneling 
spectrum is abruptly altered from a BCS-like shape to a 
much more smeared shape with substantial conductance 
below the gap energy. As dp increases further, the su- 
perconducting features in the tunneling conductance are 
strongly attenuated. Above 2.0 nm, a narrow zero-bias 
conductance peak develops which is suppressed in fields 
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FIG. 3: Superposition of five scaled conductance curves for 
(If = 0.5 to 2.5 nm. The (If = 2.5 nm curve is scaled by a 
factor of 500. 

greater than H C 2] this suggests it is related to supercon- 
ductivity, but due to its weak signal-to-noise ratio, we 
do not focus on it in the present discussion. When dp 
exceeds 3.5 nm, all recognizable features disappear and 
the normalized conductance is equal to 1 ± 10~ 4 . 

The most striking observation, though, is that between 
0.5 nm and 2.5 nm the spectra can be rescaled onto a 
single curve. That is to say, 

a(V,d F )-l=A(d F )(N(V)-l) (1) 

where N(V) is a generic, thickness independent func- 
tion and A(cIf) is a scaling coefficient defined so that 
iV(0) = 0. The N(V) curves derived from the data using 
this scaling procedure are shown in Fig. The overlap 
of the curves demonstrates the generic nature of N(V). 
A(d F ) is plotted in Fig. 0] (full circles). The straight 
line is an exponential fit with decay length d* — 0.4 nm. 
Note that for d F > 2 nm, we disregard the narrow peak 
at V — when scaling the curves. The fact that A(d F ) 
extrapolates to 1 as dp tends to zero is not required by 
the scaling procedure. It has the consequence that N(V) 
can be interpreted physically as the normalized tunnel- 
ing DOS when some minimal thickness of F material has 
been deposited below the S layer. Note also that the 
zero-bias conductance is simply 1 — A(dp) and therefore 
rises exponentially from zero as dp increases. 

An important characteristic of the generic DOS N(V) 
is that the sum rule on the DOS is satisfied within one 
percent. This not only justifies the procedure used to 
isolate the superconducting DOS, even when it is much 
weaker than the "normal" background, it also suggests 
that our tunnel barrier is not weakend by the addition 
of the CoFe layer. Such weakening would result in an 
excess current within the gap, raising the total spectral 
area above unity [24| . 
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FIG. 4: Scaling factor, A(cIf\ and normalized transition tem- 
perature, t — (T c — Min[T c ]) /T c (dF = 0), as a function of dF 
for SF bilayers. For t(dp), only the initial decay is plotted 
(dp < 2 nm). For both sets of data, the characteristic decay 
length is 0.4 nm. 



It is interesting to compare these DOS data with the 
T c (dp) data reported earlier on related samples with a 
thinner Nb layer, mentioned above. There are no evi- 
dent oscillations in A(d F ) near dp « 2 nm, where an 
oscillation is seen in T c (d F ). On the other hand, in 
Fig. we compare the initial drop in T c (open circles) 
with the scaling factor, A(d F ), derived from the DOS. 
Specifically, we plot the normalized critical temperature, 
t = (T c - Min[T c ])/T c (d F = 0). Remarkably, both sets 
of data show an exponential decrease with the same char- 
acteristic length scale. 

The absence of oscillations in our DOS data contrasts 
with a previous experimental study on a different ma- 
terial |25j. In that study, a weak ferromagnetic alloy 
(Pdi-^Nia,) was studied at two different thicknesses and 
a robust inversion of the DOS was seen. On the other 
hand, a subsequent study in which the Ni composition 
(and hence E ex ) was varied at fixed dp showed scaling 
in the DOS similar to that reported here 26]. It is nat- 
ural that these different approaches should yield similar 
results, as the relevant magnetic length, £p, is given by 
^/HD /2E ex , where D is the diffusion constant in the F 
layer. In addition, the authors of Ref. [2^| note that for 
relatively thick (dp ^> £ F ) samples such as theirs, scaling 
of the density of states with exponentially decaying os- 
cillations is predicted by the linearized Usadel equations 
appropriate to that limit. This is consistent with their 
data |2jj, but it is clearly inconsistent with ours. 

Almost any application of the Usadel formalism that 
starts with a large exchange field will predict decaying 
oscillations in the DOS. In order to prevent these oscilla- 
tions within this framework, one must include a generic 
Abrikosov-Gorkov spin- breaking parameter, Tag ^ E ex . 
Further, to obtain strictly exponential scaling of the DOS 
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(in the dp 3> £,f limit) one must take Tag to be much 
greater than E ex . In this model, d* = \JhD / 'AT ag, 
which would neccessitate Tag = 300 meV to account 
for our data. 

Although the above model gives a qualitative explana- 
tion for the two key features of our data (the exponential 
nature of A(dp) and the invariance of N(V)), it is not 
clear how accurately it corresponds to the actual physi- 
cal situation. Our thinnest sample, with dp — 0.5 nm, 
for instance, is certainly not in the semi-infinite limit, 
dp ^> d*. Thus, we would expect finite-size effects, such 
as a decrease of the gap width, which would prevent it 
from following the predicted scaling relationship. Fur- 
ther, one must wonder about the origin of Tag and the 
seeming unimportance of E ex . Previous studies on thin 
CoFe films have confirmed that near-bulk ferromagnetic 
ordering exists in films as thin as 1 nm (2^, and further, it 
is reasonable to suggest that given the constant shape of 
the DOS, none of the important physical parameters are 
changing significantly over the thicknesses we examine. 
Therefore, we conclude that a constant exchange field 
is present in all of our samples with a CoFe layer. An- 
other possible mechanism for washing out the expected 
oscillations could be lateral variations in dp j^, but we 
note that these variations would have to be vary large 
(Sdp ~ which seems unlikely. Finally, we note that 
the Usadel equation is only strictly valid in the diffusive 
limit, wherein all relevant length scales are larger than 
the mean free path. From resistivity data, we estimate 
I = 1.0 nm, which is somewhat larger than d* = 0.4 nm, 
which is comparable in size to Xp = 0.3 nm. In light 
of the totality of the considerations above, we come to 
the position that no reasonable application of the con- 
ventional Usadel theory or purely materials problem can 
account for our results. 

In summary, we have measured the tunneling DOS on 
the F side of SF bilayers and have found a sharp change in 
the shape of the conductance curves between dp = and 
dp = 0.5 nm and that for all dp > 0, the shape of the 
conductance curves is universal. The only dependence 
on dp is given by an exponential decrease in the mag- 
nitude of the superconducting signal with characteristic 
length d* = 0.4 nm. We also note a similar exponential 
decrease in T c (dp) in related samples. Finally, we note 
that we have been unable to reconcile our results with 
the conventional Usadel equation; thus, new theoretical 
considerations appear necessary. 
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